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Abstract 

We study the non-equihbrium dynamics of a system of coupled 
scalar fields in a Friedmann-Robertson- Walker (FRW) universe. We 
consider the evolution of spatially homogeneous "classical" fields and 
of their quantum fluctuations including the quantum backreaction in 
the one-loop approximation. We discuss in particular the dimensional 
regularisation of the coupled system and a special subtraction proce- 
dure in order to obtain the renormalised equations of motion and the 
renormalised energy-momentum tensor and ensure that the energy is 
well-defined and covariantly conserved. These results represent at the 
same time a theoretical analysis and a viable scheme for stable nu- 
merical simulations. As an example for an application of the general 
formalism, we present simulations for a hybrid infiationary model. 



1 Introduction 

While the renormalisation of a single scalar field in a Friedmann-Robertson- 
Walker (FRW) universe, and fermion fields coupled to such field, have been 
discussed by various authors [H El [3|, HI [5], |6] , a consistent and coherent formu- 
lation of coupled scalar field models including non-minimal gravity couplings 
and full renormalisation is still missing. As a first step in that direction. 
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we study here such a system of coupled fields in the one-loop approximation 
with special emphasis on its renormalisation in the MS scheme. This will 
allow us to obtain a consistent set of coupled equations of motion without 
divergences and numerical instabilities, which is highly suitable for numerical 
applications. Coupled scalar fields appear in cosmology in multi-field models 
of infiation and are important as well in particle physics in the context of 
grand unified theories, which contain several Higgs fields. 

The standard example of a coupled system of scalar fields in cosmology 
is of course the hybrid model of infiation which has been introduced by 
Linde [3 El |9] and whose coupled dynamics at the end of the infiationary 
phase has received wide attention [T0l[IIl[T2l[T3l[Ta|15l|16lll71|18l|l^ 
EDI ED [221 ESI [21I2S] (for a recent review see P]). Unlike other models 
for infiation and for preheating, where one has just one "infiaton" field with 
vacuum quantum numbers coupled to various other boson and fermion fields, 
this type of models has two scalar fields with can acquire a time dependent 
vacuum expectation value. So these fields and their quantum fiuctuations 
can mix dynamically. 

The renormalisation of such a coupled system of bosons out of equilib- 
rium has been considered previously by Cormier et al. [22] in the one-loop 
approximation, and in Ref. j27j in the Hartree approximation. The one-loop 
approximation includes the quantum back-reaction of the quantum fiuctua- 
tions on the classical fields, the Hartree approximation includes in addition 
the backreaction on the quantum fiuctuation themselves. This approxima- 
tion has been considered as well by Bastero-Gil et al. [13] with a momentum 
cutoff for the quantum fluctuations. These publications do not include the 
coupling of the flelds to gravity, which may be a reasonable approximation 
in the period of preheating after inflation. However, when describing the 
evolution during inflation and the end of inflation, the inclusion of the cou- 
pling to gravity could be important and it is necessary for consistency, since 
those couplings are generated by quantum corrections. Moreover in recent 
years, such non-minimal gravitational coupling has drawn a lot of attention 
in the context of realising inflation within the Standard Model [28] 129] and 
it is surely an interesting issue to include it also in the case of many flelds. 
In those models it became clear pretty soon that quantum effects on the 
inflation potential cannot be neglected [301 ED |32] and it may be the same 
also for the hybrid case. 

The one-loop renormalisation for the out-of-equilibrium evolution of a 
single scalar fleld in a FRW universe has been considered by several au- 
thors [D El El ESI EH ESI ES] • The technical approach which we will use here 
is based on Ref. [36], which itself uses a formalism developed in Ref. [37]. As 
in these references we establish here a " renormalised computation scheme", 
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i.e., the renormalisation is taken into account with the same rigour as in more 
formal approaches, and at the same time the formahsm can be implemented 
efficiently into a numerical code. In particular the renormalisation is inde- 
pendent of the initial conditions, though the divergent momentum integrals 
depend on the initial effective masses of the system. The energy-momentum 
tensor is in this formulation always finite and covariantly conserved. 

The paper is organised as follows: in Sec. [2] we recall some basic relations 
of FRW cosmology, with special attention to an extension beyond space-time 
dimension n = 4 and with the inclusion of higher dimensional gravity coun- 
terterms, which are needed in order to properly renormalise the coupled fields 
system; in Sec. [3] we define the general model whose dynamical equations 
we will formulate in the one-loop approximation; in Sec. IHwe discuss in de- 
tail the regularisation and renormalisation of the equation of motions and of 
the energy-momentum tensor. Finally in Sec. [5] we give few selected results 
to show the applicability of this scheme and the stability of the numerical 
results and then conclude. More exhaustive numerical results will be given 
in a forthcoming publication. Various technical details can be found in the 
Appendices RllUl 



2 Short review of FRW cosmology 

We consider here a spatially isotropic and flat space-time with n — 1 spacial 
dimensions, which can be described by a FRW metric with k = curvature 
parameter as 

ds^ = dt^ -a^{t)dx^ , (2.1) 

where dx^ = Yl^=o(^^i- The effective action of the coupled gravity-scalar 
fields system is given by 



2 

+ Ag.^ + 6aR^ + 6 PR "^R + 6jR ""^^^R , 



+ J d^'x Ci<^,,g^,,R) . (2.2) 

The first two terms are the Einstein-Hilbert action and a cosmological con- 
stant term, which allows to tune the energy density of the ground state to 
zero. We have introduced also all the gravitational terms up to dimension 
four for later convenience. The last term is instead given by the action for 
the matter fields, in our case a system of coupled bosons described in the 
next section. 
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The time evolution of the a{t) is governed by the Einstein's field equation, 
which for the Einstein-Hilbert case 5a = 5/3 = ^7 = is given by 



G^, + Ag^, = -8nG <T^,,> . (2.3) 
where the general expression for energy-momentum tensor is 

T,Ax) = -^^-^p^ . (2.4) 

In quantum field theory renormalisation requires to introduce a countert- 
erm action. Besides the usual mass and coupling constant counterterms the 
coupling to gravity induces divergences which require the inclusions of the 
higher curvature terms written in Eq. (12.21) and therefore a modification of 
the original Einstein field equations. These then take the form 

{l + 6Z)G^, + 5a^'^H^, + 6p^^^H^, + 6^H^, + {A + 6A)g^, = -8nG <T^,> . 

(2.5) 

The tensors ^^^H^,y,^'^^ H^y and Hf^^ are related to the variation of the higher 
curvature terms B?^ R"^ Rap and R"^'^^ Rap-yS- Since gravity is not quantised, 
all divergences arise, at least technically, from the quantum fiuctuations of 
the matter fields in the energy-momentum tensor. Therefore one may prefer 
to include the counterterms into a redefinition of the right hand side of the 
Einstein equations via 

T^u ^ T^u + ^ZG^y + da^^^H^y + 6~P^^^H^, + 6^H^, + (A + 6K)g^y . (2.6) 

Here 5Z = 6Z/8ttG and similarly for the other terms. 

For the various tensors we use the conventions of Ref . |38j : the Einstein 
curvature tensor is defined as 

G^u = Rfiu ~ -g^luR', (2-7) 

the Ricci tensor and the Ricci scalar are 

Rfif = R^vx^ R = g'^'^Rfj.u ! (2-8) 
where the curvature tensor is given by 

= '^■y^ap ~ daT^p + r"^r^^ — r^^r^^ . (2.9) 
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The higher curvature tensors are defined as 
1 6 



-9 Sg 



2R. 



1 



1 6 



2g UR — -g ij^y R + 2RR ^i, , 



(2.10) 



2R 



-9 



j (Tx^/^R'^^R 



a/3 



UR 



R -fiu — -g fiv ni? — UR ^y — -g ^^R °'^R + 2R "'^R a^ipv 



H 



-g Sg 



j (Tx^^R'^^^^R 



(2.11) 



fj,aj3u a/37 



4aR^, + 2R 



For the case n = 4 in conformally flat space-time one has 



fJ,U 



(2.12) 



(2.13) 



The exphcit expressions for all these tensors in FRW space-time are of 
course well-known. However, we will use dimensional regularisation, which 
implies that we have to perform all the computations for general n = 4 — e. 
We therefore recall the explicit formulae for the relevant tensors in flat n 
dimensional FRW geometry in Appendix |X1 

For the spatially isotropic FRW universe the Einstein fleld equations re- 
duce for the time-time component and trace to the Friedmann equations 



Gtt 



-SvrGp 



-87rGT; = -87rG(p-(n-l)p) 



where p is the energy density, and the pressure p is given by 

p=(p-T;)/(n-l). 
. The covariant conservation of energy takes then the form 

'^^={n-l)H{p + p) = H{np-Tli) . 



(2.14) 
(2.15) 



(2.16) 



(2.17) 
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3 Coupled fields in curved space-time 

We here consider the quantum field theory of N coupled scalar fields. More 
specifically we consider a class of models whose Lagrangian is given by 



with the potential 



N N N 

vm = ^ E ^"^^ + i E E ^^.-^^'^1 • (3-2) 

1=1 i=l j=l 

Instead W{R, 4>) describes the coupling of the scalar fields to the curvature 
scalar: 

W{R, $) = J] |i? . (3.3) 

i 

The are dimensionless coupling constants that have no counterpart in the 
fiat space theory. As we will see later there is an obviously preferred choice 
= 1/6, called conformal coupling. The value = is called instead 
minimal coupling. 

3.1 The equation of motion in the one-loop approxi- 
mation 

We split the fields in two parts, the expectation values and the quantum 
fiuctuations around it 

^i = cPiit)+^^{t,x) , (3.4) 

where the classical field (pi(t) is assumed to be homogeneous in space. The 
classical part of the Lagrangian density retains the form 

-|E^^^'^M-^| ■ (3-5) 

ij ) 

The first order in the fiuctuations vanishes (the expectation value of the 
fiuctuations is zero). The fiuctuation Lagrangian density of the second order 
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in fluctuations 



\ i i 



(3.6) 



Now we can write the equations of motion for classical fields: 

4>i + {n-l)H<j)i + {m\ + i,K)(\), 

+ 5^ < V'l >)'/'i + 2 < i^ii^i > <Pj] = , (3.7) 

j 

and for the quantum fluctuations in the one-loop approximation: 



The latter equation can be written in the form 

j 

where the mass squared matrix is given by 



or, more explicitly, 

Ml 



(3.9) 



(3.10) 



(3.11) 
(3.12) 



The expectation values of the quantum fluctuations can be expressed in terms 
of equal-time Green's functions via 



< ViV'i >= -iGij{t,x;t,x) , 
where the Green's functions satisfy 

, 1 



(3.13) 



^S,,^S^-\x-x') 5it-t'). 
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(3.14) 



We have introduced here the set of equations for the classical fields and 
for the fluctuations in a heuristic way. A rigorous derivation can be found 
in Ref. [39], based on the Schwinger-Keldysh or closed-time path formalism 

gniii]. 

The differential operator in Eq. fl3.14p is non-Hermitian due to the term 
{n — l)Hd/ dt, and so is the matrix valued Green' s function G. This problem 
can be solved by introducing conformal time 



T 



* 1 
dt' 







a{t') 



(3.15) 



We can then introduce the dimensionless rescaled fields and their fluctuations 



T 



-n/2+1 



as well as the rescaled Green's function 

G,,(r,x;r',x') = a'^/^-^t) a^'^-\t') Gi,{t,^;t' ,^') 



(3.16) 



As a function of conformal time, the Hubble parameter, the Ricci scalar and 
the derivatives of the fields are given by 



H 
R 
f 
f 



2{n - 1)— + {n- l){n - A)H^ 



-n/2 



f- 



n 



-aHf 



(3.17) 



-n/2-l 



o n — 2 



—a 



4(n- i; 



-!)(»- 2) 



Rf 



where / stands for any field and the prime denotes the derivative with respect 
to conformal time. 

Then the equations of motion of the classical fields become 



where 



n-2 
4(n - 1) 
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= 0, 
(3.18) 

(3.19) 



We see here that for n = 4 and conformal couphng = 1/6, the curvature 
term disappears from the equation of motion of the fields, which reduces to 
the form it has in flat Minkowski space if also = 0. 
The new two-point functions satisfy now 



(5(r-r')(5"-^(x-x') 



with the effective masses of the rescaled fluctuation fields 



Miir) 



The equations for the fluctuation fields in conformal time become 



V^/' + V2^, + 5^^^^.(r)^, =0. 



(3.20) 



H + - ^n)R) a' + a'-" ( SXj' + J2 K^] ) '(3-21 



(3.22) 



(3.23) 



We expand the fluctuation fields in terms of mode functions /"(t; k) 



(27r)' 



which satisfy the mode equations 

/r(r; k) + k^ftir- k) + -M^,.(r)/;(r; fc) = . 



(3.24) 



(3.25) 



The latin subscripts refer to the field components, i = 1, . . . , and the greek 
superscripts a = 1, . . . , refer to the independent solutions of the system 
of coupled differential equations 03.251) . We will specify them below, by a 
suitable set of initial conditions. 

The Green's functions can be expressed by their Fourier components 



A_le— ')G,.(r,r';fc). (3.26) 



These Fourier components can be rewritten in terms of mode functions. They 
read 



G,,(r,r';fc) = J^W^^ 

a/3 



fnr,k)ff{T',k)eir-r') 



+fnr',k)ff{T,k)e{r'-r) 



(3.27) 
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where Wajs is the Wronskian matrix of the system of solutions: 



(3.28) 



Therefore the expectation values of the fluctuation fields are given by the 
fluctuation integrals: 

/^n — 1 / 
77r-^,J2^i^^e{ftft)- (3.29) 

^^^^ a/3 

The Wronskian matrix is determined by the initial conditions of the system 
of fluctuations which will be discussed in the next subsection. 



3.2 Initial conditions 

The choice of initial conditions for the quantum system is very important 
for the renormalisation of the equations of motion. In order to satisfy the 
canonical commutation relations for the creation and annihilation operators 
and choose a Fock space basis, we first diagonahse the mass matrix Mjj at 
the initial time. The eigenvalues rri^Q, a — 1,N then define N independent 
free particle Fock spaces whose product we choose as our initial state. The 
eigenvalues and eigenvectors of the initial mass matrix A4fj(0) are defined by 
the equations 

det{A^|(0)-my,,}=0 (3.30) 

and 

-<o/ro + ^',(0)/;o = 0. (3.31) 
The eigenvectors are orthogonal, we choose them real and normalised as 

E /"o/fo -^ap, E fS>fZo = ^im . (3.32) 

I a 

We now specify the mode functions introduced in the previous section by the 
initial conditions 

fnO,k) = fro, (3.33) 
fr'{0,k) = -tn^ofto^ (3-34) 

with the frequencies ^lao{0) = a/Zc^ + ^qo- The Wronskian matrix of the 
mode functions at time r = is given by 
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It is time-independent by construction. With these initial conditions the 
fluctuation integrals are given by 



///"-It- 1 

These initial conditions correspond to a quantum state which is denoted as 
the adiabatic vacuum. For such an initial state the energy-momentum tensor 
is found to be singular at initial time |12]. This problem can be solved by a 
suitable Bogoliubov transformations of the adiabatic vacuum. For a coupled 
system this Bogoliubov transformation was derived in Ref. |13]. We recall 
the essential features of this approach in Appendix [Bj As discussed in Ref. 
|12] this modification of the initial state does not affect the renormalisation, 
which is our main subject here; so we continue to work with the adiabatic 
vacuum as the initial state. 



3.3 The energy- momentum tensor 



In order to write the Einstein's equations we need the expression of the 
energy-momentum tensor for our system of coupled field in Eq. (13. ip . Gen- 
eralising straightforwardly the case of a single scalar field given in Ref. [38], 
we have 



T, 



fJ,U 



(3.37) 



In the Friedmann equations only Tu and appear and their classical parts 
are given by 



T'C 



rjicl fl 



with 



E 

i 

E 

i 



l-- + 2(n-l)e. 



• 2 



+ 2(ri-l)ei 4>i + {n-l)H(i)i 



(3.38) 

(3.39) 
(3.40) 
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Using conformal time and the rescaled fields we find 

^ 4 



2 2 12\ 



n - 2 



j4— n 



(3.41) 



For the fluctuation contribution to the energy density we obtain instead 



X 



rE 



1 



2fi 



a" 7 (27r)' 



(3.42) 



+(r!-l)5^(6-f„) 



72 — 2 9 r, i?a^ 



2 n - 1 



\fi 



a\2 



The last term, proportional to which has no analogy in the classical 
energy-momentum tensor, is not genuine but compensates an analogous term 
in M-fj- It is convenient to rewrite T^^ as a function of the fluctuation integral 

as 



1 



2a" 



n — I \ ^ ,^ 



aH—J'ii 

CLT 



n — 2 9 Ra^ 
——a^H^ + 



2in - 11 



where we have introduced the "potential'' 



V,,{r) = MUr) - Mm 



(3.43) 

F 

•> II 

(3.44) 



and where 



t: 



iq,kin 



1 1 



a" J (27r)' 



(3.45) 



^Though we have introduced conformal time we continue to use the time-time compo- 
nent of the energy-momentum tensor in standard time as made explicit by using time 
indices t instead of or r. 
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is the kinetic or free-field part of the energy density. 

The classical part of the trace of the energy-momentum tensor becomes 



(3.46) 



n 



For the quantum contribution one finds instead 



a" J {2rr) 



d''~^k ^ 1 



2n 



X 



n-2 „ d 



2 dr ' 



n 



aH-\m' + -{{n-2YH' -R\\f, 



n — 1 



COL\2 



(3.47) 



Using the equation of motion for the mode functions, Eq. fl3.25p . the first 
three terms in the bracket can we rewritten as 



i/r? - k'\ft? - Myrff = i/r? + ffft" = ^^i/rP • (3-48) 

We therefore may express the entire quantum contribution to in terms of 
the fluctuations integral J-'ij as 



n — 1 



ar"^ dr 



(3.49) 



+ ^ ( {n-2fH^ 



n 



n — \ 



R 



4 Renormalisation 

In the previous sections we have introduced the fluctuation integrals and 
expressed the equations of motion and the one-loop contributions to the 
energy-momentum tensor in term of these quantities. These definitions have 
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only formal character, as they involve divergent momentum integrals. As in 
the case of renormalisation in Minkowski, we have to introduce regulators 
in order to render the integrals well-defined, and counterterms to absorb the 
divergences when the regulators are removed. We will use here dimensional 
regular isat ion with the space dimension 3— e. Furthermore, it is convenient to 
separate the limit e — )■ entirely from the numerical computations by adding 
and subtracting appropriate analytical terms, which reproduce the divergent 
behaviour and render the numerical integrations finite. Such a scheme was 
set up in Ref. (3^, and we will use it here. We recall its main features and 
the formulae needed here in Appendices O, O and [El 

Before we discuss renormahsation we have to discuss a subtle issue (see 
also Ref. We have written all equations for general n = 4 — e, and the 

equations of motion ensure that the energy-momentum tensor is conserved 
in n dimensions. Ultimately we have to take the limit n — )• 4, or e — )■ 
and this will generate terms that behave as 1/ (n — 4) from the divergent 
integrals and consequently in the counterterm action. On the other hand the 
equations of motion and the energy-momentum tensor are equal to their 4- 
dimensional form only up to terms of order (n — 4) . These terms will survive 
the limit n — )■ 4 whenever they are multiplied by a factor l/(n — 4). This will 
result in "additional finite terms", which have no counterpart in quantum 
field theory in Minkowski space. So our prescription is: we first formulate 
the full renormalised theory in n dimensions and take the limit n — > 4 only 
at the end. Another prescription one may think to use, is consider n ^ A 
only in the divergent integrals, while keeping the equation of motion and the 
energy-momentum tensor in 4 dimensions, but in that case we cannot rely 
on energy-momentum covariant conservation during the whole regularisation 
procedure. 



4.1 Renormalisation of the equation of motion 

The equation of motion contains logarithmic and quadratic divergences due 
to the fluctuation integrals. Using the analysis of Appendix O and in partic- 
ular Eq. flC.12p . their divergent behaviour can be understood and it coincides 
with that of the divergent integral 



-pat 



■div 



d'^'^k ^ 1 



(2vr) 



fa ra 
JioJ jO 



V, 



(4.1) 



where the potential Vap has been introduced in Eq. (IC.lOp Here and in the 
following we use a slight change in notation, without introducing a new sym- 
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bol: we incorporate the prefactor which appears in front of the divergent 
integrals into the ra- dimensional integration measure. 

Using the formulae for the various integrals given in Appendix [El one can 
obtain analytically the dimensionally regulated form of the previous expres- 
sion as 



where = |— 7+ln47r. In order to cancel this divergence in the equations of 
motion we have to add an appropriate counterterm to the Lagrange density. 
In the case of Minkowski space-time, it was found in Ref. ^21 that the 
counterterm Lagrangian 

^mnk. = -SCMlMk (4.3) 

with 



takes account of all divergences. We work here in the MS scheme and sub- 
tract only the divergent part of the fluctuation integral. Note that for con- 
stant masses and fields, this is a contribution to the renormalisation of the 
cosmological constant of order mf and also contains mass and coupling 
renormalisation terms due to the interaction part of the effective mass matrix 
Ml. 

Taking into account the factor y/—g = a"" (using conformal time), and 
noting that the Aijj differ by a factor from Mfj, the equivalent Lagrangian 
for general FRW is 

^ -a-'^5(^MlMl . (4.5) 

When discussing the energy- momentum tensor it will be convenient to have 
this Lagrangian decomposed into the usual mass, coupling constant and other 
gravity-related counterterms. This is presented in Appendix [Fl 
Introducing this counterterm Lagrangian and given that 

d(j)i V 

= ASCJ2[^MI\u4>i + XhM1i^ , (4.6) 
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the equations of motion for the classical fields become 
+ (m^ + (e. - \)R)a% 



(4.7) 



. 



It follows from Eq. (14.21) and the renormalisation procedure that the combi- 
nation A5C,J^'jj + J^ij is finite. This is entirely analogous to the case without 
coupling to gravity. There is, however, one important difference: in addition 
to the renormalisation scale the scale factor a(r) appears in the prefactor 
and this modifies the finite terms. In order to be able to identify the different 
contributions to J^ij separately we write 



add 



where the finite part of the fluctuation integral is given by 



■fin 



TTSub I -p{t 
•'ij T" ''ij ■ 



(4.8) 
(4.9) 



Here we define 



■sub 



27r2 



^21],o 



V„«(r) 



(4.10) 

as the subtracted fluctuation integral, that is finite, but in general has to be 
computed numerically. From J^iJ^, Eq. (14.21) . we have 



E 



m 



aO ra ra 



In 



- 1 



+ 5^ fiofjO 



167r2 



/30 



In 



■^0 



- 1 



m 



aO 



(4.11) 



as finite terms that are left over after removing the divergent part. However, 
among those we have identified separately 

1 



X-add _ 2 p r 

which arises from the 0{e) term in the expansion of Aiff- 



Ml = {m\ + (^i - l/6)i?) a^ + a^— SA^ 



r 2 




(4.12) 



+ -a2i? + 0(e2) 
do 



(4.13) 
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This is the only term that has no usual AD counterpart and is generated 
from taking the n-dimensional FRW model. It gives a finite contribution to 
the scalar field coupling to gravity C,i and it may be in principle absorbed 
into a more general renormalisation scheme than MS. Note that such term 
survives even in the case of 4D conformal coupling = 1/6. 

Finally, the renormalised equation of motion for the classical fields reads 

4>■+{m^ + {^^-hR)a'4>^ (4.14) 



potential Vq,^ in J^fj"' is computed in 4 dimensions. Likewise, the fiuctuations 

are computed using the 4-dimensional reduction of Ai'^j in their equation of 
motion. 

Note that the occurrence of In a(r) in the finite terms A4fj is not nec- 
essarily a small correction. Neither is the "potential" Va/s small: the terms 
{mj + — 1/6)R) appear multiplied by a^. Indeed the qualitative role of the 
In a(r) term in J^/j is to compensate a logarithmic increase of the subtracted 
integral J-'ff^. This is illustrated by an example in Appendix iGl 

4.2 Renormalisation of the energy-momentum tensor 

We can now proceed in the same way to the renormalisation of the energy- 
momentum tensor. In Sec. 13.31 we have separated the quantum part of 
the tt component of the energy-momentum tensor into a "kinetic" part and 
some further contributions which can be written in terms of the fiuctuation 
integrals. The divergences of the latter have been discussed in the previous 
section, so it remains here to analyse the kinetic part. Using the results of 
Appendix O and in particular Eq. (jC.181) we obtain 



a" J (27r)"-i ^ 2VL^Q 



1 
2' 



02 I _|/,a'l2 



(4.15) 



The divergent behaviour of the first term in the bracket is obvious, the one 
of the second term is given by Eq. (]C.20p . The third term originally appears 
with the mode functions /f ; the present form, with the reduce mode functions 
is obtained by using that the parenthesis (-Mfj(O)— m^Q^jj) vanishes when 
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multiplied with or /"q due to the eigenvalue equation fl3.3ip . Still the 
term could lead to a divergence arising from the second order contribution 
can convince oneself that this leading term vanishes. The 



expansion of hi and contains factors /^^ or /^q (see Eqs. fIC.Sp and 

fIC.Qp ): the parenthesis then reduces to ('^^o ~'"^ao)^u ^^^^ found to be 
multiplied with a term that is symmetric in a and /3. So, using the integrals 
given in Appendix [El we can write 



-iq,kin 



1 



647r2a^ 



a aj3 



q, kin, fin 



(4.16) 



Here 

rT-iq,kin,fin mq, kin, sub , rT-iq,kin,ft 

^tt — -'-tt -'-tt 

where the subtracted integral is given by 
1 f k'^dk 



q, kin, sub 



tt 



i_ rk^^j_ i|,„,|2_v- 

+ ^(-M',(0)-my,,)Re(/.r/^r) 



(4.17) 



/30j 



(4.18) 



and the finite terms of the regularised integrals are 




1 "^aO 



In 



/30 



2 2 2 



- 1 



(4.19) 



In both of these equations it is understood that n = 4. 

Besides T/^''^™ the tt component of the energy-momentum tensor contains 
further quantum contributions, which have been expressed by fluctuation 
integrals, see Eq. fl3.43p . Collecting all these finite pieces we define 



t: 



iq,fin 



pq, kin, fin 
-tt 



■fin 



—fin ( „2 zt2 I \ jrfin 



(4.20) 



where we have already set n = 4. 
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For the singular terms we find instead 



IGvr^ dr 

Here it is understood that everything is still defined for general n. The 
counterterms are (see Eqs. fl2.6l) and fl3.4ip ) 



^ i 

+2(n - 1) 5^ '^^i - "^^a^H^^j 0^ 

+^ ^ 5A,,0,20^2 I + + 6ZGu + ^a^'^i^it • (4.22) 



Adding the singular part and the counterterms one finds that the divergences 
are cancelled, but, as in the Cc 
in the limit n — > 4. These are 



are cancelled, but, as in the case of J-ij some additional finite terms remain 



n^4^ " 967r2a' 



2 



(4.23) 



It is understood that n is set equal to 4 in yVl^j and ^^^Hu- The result agrees 
with the one for the single- field case found in Ref. [36]. There T^'^'^'^ was 
in addition presented in an expanded form, using the explicit expressions 
for Ad'^i and ^^^Htt. Here we refrain from displaying such a rather lengthy 
formula. We should mention that the calculations for checking the cancel- 
lation of divergences and for obtaining T^''^'^'^ are quite cumbersome, they 
have been performed using the computer algebra code REDUCE |44]. Note 
that these additional terms are time- dependent and similar time-dependent 
pieces are contained in the subtracted integrals eq. fl4.18p . as it happens for 
the In a(r) terms in J^ij. Both together ensure that the energy-momentum 
tensor is covariantly conserved. 



19 



There is still a further finite contribution to the energy-momentum tensor, 
the conformal anomaly, which can in any case not be renormalised away [5S] . 
As discussed there and in Ref. [36] it cannot be determined within our 
framework but has to be taken over from a more general analysis. It is 
obtained by setting, for N scalar fields, 

5/3 = -5^ = ^ . (4.24) 
^ ' 28807r2 ' ^ ' 



As 



H^y in 4 dimensions, the n 



4 part of these tensors cancels, and 
so does the singularity. But there is a finite remainder from the extension to 
n 7^ 4. We obtain 



T-iano 



lim - 

n->4 28807r2 



N 



28807r2 



R' 



H— + RH' + —R' 
a 12 



(4.25) 



Notice that these additional terms are vanishing for the case of pure de Sitter 
or radiation dominated expansion, so they are negligible in the infiationary 
phase, but they could have an effect in the reheating phase. 

Finally, the renormalised tt component of the energy-momentum tensor 
is given by 

+ T^'l'^'' + T.'l''"^'^ + ° . (4.26) 

\'S.'S\ ^, the unrenormalised quantum contribu- 
tion to the trace of the energy-momentum tensor can entirely be expressed in 
terms of the fluctuation integrals J^ij and derivatives thereof, see Eq. fl3.50p . 
So the flnite part is simply given by 



T-iren 

-'-tt ~ 

As we have seen in Sec. 



T-iq,fin /i 



+ 2H' 



1. 
6' 



■fin 



2aif-^J-f 

dT 



R\T\ 



fin 



+ 



fin 



(4.27) 



Likewise, the divergences of ^ are related in a straightforward way to those 
of the fluctuation integrals. Therefore, the part that is singular as — 4, is 
given by 



7^q,sing ti 



167r2a" 



\n — 



^\([n- 2fH' - -^R 1 
n — 1 



+ Y,^,Ml\ ■ (4.28) 
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We have to add the appropriate counterterms which have the form (see Eqs. 

and dsan])) 



1 



b[--{n-2)aH<P, 



(4.29) 



+n 



i 

i 

i ij 

+6ZG'', + s&'^^^m . 

We again find, using REDUCE, that the divergences cancel in the sum of 
j'q-singAt rpq,ctij.^ There remain finite terms, however, as for T^. These are 
given by 



'^■^^^-2ai7^ + 2a2ij2_;^2 



dr 



(4.30) 



12^^ + l(i?-18i/2)^2 -^2_l(i)^,. 
;3 (ir ^ " 18 36 ^ 



a-" 



These agree, for the single-field case, i.e., = 1, with Eq. (7.28) of Ref. 
[36] . except for a misprint there: {R + 18if^) should read {R — 18H^), as 
here. The expanded form (7.29) in Ref. [36] is correct. We again refrain 
from presenting such an expanded form. 
The anomalous part is obtained as before 



A* 



N 



lim -L, ((2)ff/;- ff/;) 



w4 28807r2 
fR" 



, , , + 2/7— + 2/2i72 _ i2H^ 
28807r2 V a2 a 



(4.31) 



Finally, we have to add the classical part, Eq. fl3.39p . with n = A. So we 
finally have 



Tiren fi ^cl fi | _i_ ^q,fin /i , ^q,add fi , rpa,no /i 
II —-I'll n=4 "T // T -t ,, ~\~ 1 ., 



(4.32) 
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Figure 1: The evolution of the scale parameter a{t) for set I. 

5 Simulations: the hybrid model 

While we have discussed here the renormalisation of a relatively general 
model of N coupled fields, the original incentive of our investigation was 
the application to the hybrid model of inflation. In that case we have a 
system of only two fields and we can write the Lagrangian as 

^ = V^{l9'''d,<j>d,<j>+^-g^''d,xd.X-lm'<p'-^a{x'-v'f 

-^A0V-^i?(0^ + X^)} , (5.1) 

Here v is the vacuum expectation value of x- Furthermore, we have to 
make use of a bare cosmological constant A = av"^ so that the minimum 
of the classical potential has the value zero. Even when choosing the MS 
prescription for the renormalisation of the masses and couplings we have to 
make sure that the zero point of the quantum fluctuations at the minimum 
vanishes, so we have to add a suitable finite counterterm 6 A. Of course this is 
"fine tuning" and represents the practical aspect of the cosmological constant 
problem. We have restricted our simulations to the case G = ^2 = 1/6, the 
conformal couplings. This considerably simplifies the dynamics, and this 
choice would be natural in a supergravity theory. 
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Figure 2: The evolution of the inflation field (f){t) for set I. 

As parameter set I we choose (in Planck units) v = 1, m = 0.01, A = 10~^, 
a = 10~^. The initial values of the fields are 0(0) = 4 and x(0) = 0.01. In 
Fig. [T]we show the evolution of the scale parameter a{t). We see an initial 
exponential expansion, followed by a power behaviour a ~ t"'®^ as expected 
for a matter dominated universe. The exponential behaviour is not due to a 
slow roll of the "inflaton" field which oscillates with a rapidly decreasing 
amplitude, see Fig. |5J Rather it is the field displayed in Fig. [3] that 
takes some time before falling from the metastable maximum to the stable 
minimum. Quantum fluctuations build up after the field x starts to fall from 
X = Otox = v = l. Of course both the classical as well as the quantum 
energy density decrease with the expansion of the universe. We therefore 
plot in Fig. mthe ratio p^/p of the fluctuation energy and the total energy 
densities. The fact that the classical energy density is not converted entirely 
into fluctuation energy is typical for the one-loop approximation in which the 
classical fields and the fluctuations remain coherently coupled. 

The covariant energy conservation, Eq. (12.171) with n = 4, is fulfilled with 
a relative error less than 1 ppm. 

We also display, in Figs. O and O the integrands of the fluctuation in- 
tegrals J^ii and J-22, including the factor /c^/(27r^). For J^n we display the 
unsubtracted integrand, the integrand with the zeroth order term removed, 
and fully subtracted with zeroth and first order in V removed, see Eq. (14.101) . 
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Figure 3: The evolution of the waterfall field x{t) for set I. 
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Figure 5: The integrand of the fluctuation integral J^u: (a) without subtrac- 
tion; (b) after removing the leading term; and (c) with zeroth and flrst order 
in V subtracted, respectively, for set I; crosses: the numerical results; solid 
lines: in parts a and b these indicate the asymptotic behaviour p and j9~^, 
respectively; in part c the solid line reproduces the approximation Eq. (IG.Sp . 
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Figure 6: The integrand of the fluctuation integral T22 for set I at t = 6306; 
the crosses are the numerical results, while the solid line is the asymptotic 
approximation Eq. (IG.Sp . 

The subtraction is actually performed using the functions /i" introduced in 
Appendix [O in order to avoid small numerical differences. At large momenta 
the integrands are seen to behave as p, and p"^ respectively. For the fully 
subtracted integrand we compare the numerical result with the approximate 
formula f lG.Sp . The agreement is almost perfect: the time at which this inte- 
grand is extracted is t = 3154. There the field has reached its asymptotic 
value already and corrections due to the field x are small. The integrand for 
J-22 is plotted at t = 6306. There the field x has come close to its asymptotic 
value 1, but is still oscillating around it. Here the approximate formula (IG.SP 
describes well the asymptotic behaviour, but one notices a strong peak at low 
momenta p ^ .1. This is in agreement with the fact that a significant part 
of the classical energy has been converted into fluctuation energy, as already 
displayed in Fig. HI Obviously it is the waterfall field whose fluctuations at 
low momenta have been strongly excited. 

6 Conclusions 

We have presented here a general analysis of the renormalisation of a set 
of coupled scalar fields in a spatially flat FRW universe with non-minimal 



26 



gravitational coupling and higher curvature gravity terms in the one-loop 
approximation. By a suitable subtraction procedure we obtain the equations 
of motion for the classical fields and the energy-momentum tensor in renor- 
malised form. The latter appears of course on the right-hand side of the 
Friedmann equations. The counterterms are found in the standard form of 
quantum field theory, they do not depend on numerical cutoffs nor on the 
initial conditions. The expressions to be evaluated on the computer are fi- 
nite from the outset, in particular all the subtracted integrals are finite and 
numerically well-behaved. 

The analysis of the divergences used here differs from the procedure used 
often in the context of quantum field theory in curved space: the adiabatic 
subtraction, using the eikonal expansion. We have not seen this approach be- 
ing applied to a coupled channel system; as far as we see this meets problems 
of principle: while ordinary functions and their derivatives commute, this 
does not hold for matrix valued functions. Our renormalisation procedure 
does not suffer from these shortcoming and can be applied to the general 
case of coupled scalar fields. We expect it to be easily generalised to contain 
fermionic fields as well and be viable also for supersymmetric models. 

The example of coupled scalar fields that has attracted most attention is 
the hybrid model of inflation. We have now implemented this renormalisation 
procedure for this specific case and obtained first encouraging results: the 
energy-momentum conservation is recovered in our scheme to high precision 
and our asymptotic formulae for the divergences describe the divergent fluc- 
tuation integrals with very good accuracy. In the renormalisation procedure 
we obtain additional finite terms, some of which arise from the n dimensional 
FRW setting and have no 4D counterpart. These terms are in some cases 
large and time-dependent, containing ln(a(r)), but they cancel correspond- 
ing time-dependent terms in the fluctuation integrals and they are therefore 
crucial for obtaining consistent results, as is discussed in Appendix [G] 

We have shown for a very simple case that our renormalisation proce- 
dure is solid and gives very stable numerical results. We plan to extend our 
study to the inflationary phase and to different hybrid models in a future 
publication. 
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A Tensors in n-dimensional flat FRW 

geometry 

The extension of the FRW geometry to n dimensions is done in a specific 
way: the three-dimensional space of FRW geometry is conformal to three- 
dimensional Minkowski space; it is this three-dimensional space that we ex- 
tend to3 — e = n — 1 dimensions in the same way as it is done in the usual 
dimensional continuation. 

For the Christoffel symbols one finds 



r*,-aa<^., ; ^i, = ^4 = --<5^ (A.i) 

(Jj 

The non-vanishing components of the Riemann tensor are 

R %k = aa Sik , R %t = -aa Sij , (A.2) 



Kk-^Si , Rljt---J, (A.3) 
for those of the Ricci tensor one finds 

Rtt = {n-1)- , Rij= [-da - (n - 6ij . (A.4) 
This leads to the Ricci curvature scalar 



R^2{n-l)- + {n-l){n-2)(-) . (A.5) 

Expressed in terms of Hubble 's constant 

H{t) - ^ (A.6) 
a{t) 

it takes the form 

R^{n-1) (2H + nH^^ . (A.7) 
The relevant components of the Einstein tensor become 

Gu = J^1^^}^H^, (A.8) 

G', = -^R- (A.9) 
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The time-time components and the trace of the tensors '•"^-ff^j/ are given 

by 

^^^Htt = -GHR + ^R^ -6H^R (A.IO) 
+ {n - 4) (^-2HR -{n + 1)RH^^ , 

^"^^Hu = -2HR+-R'^ -2H^R+(n-4) (--HR (A.ll) 

6 \ 2 

--{n + 2)H^R+-{n-l){n-2fH^ 



24(n - 1) 4 



Htt = -2HR+^R^ -2H^R (A.12) 



+ (n-4) 



WH'/^ = -6R-18HR (A.13) 
+ {n - 4) 1^-2/2 - 2{n + 2)HR - ^i?^ 

(2)^/^ = _2i?-6i7f?+ (n- 4) 1^-^/2- ^(n + 3)i/i? (A.14) 

+-{n- 2fH^R - -n{n - l)(n - 2fH^ 



8(r2-l) 4 

H^^ = -2R-6HR+(n-A) {-2HR- (A.15) 

^ \ 2[n — 1) 

+ {n - 2)H^R - ^n(n - l)(n - 2)i7^ 

B Removing the initial singularity 

The initial conditions we have discussed in Sec. 13.21 lead after renormalisa- 
tion to singularities in the fluctuation integrals when r — )■ 0. They contain 
terms proportional to Va/3(0) and Va/3(0) which behave as rln[(mQ,o + "^/3o)''"] 
and r^ln[(mQ,o + 'm^Q)T\ respectively. This does not affect the equations of 
motion, but the energy-momentum tensor, which contains the first and the 
second derivatives of the fluctuation integrals, that are infinite as time goes to 
zero. This means that these initial singularities will appear in the Friedmann 
equations as well. In Ref. it was shown how to remove such singularities 
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by a modification of tlie initial state. Here we will just briefly repeat these 
results. 

The general concept of the Bogoliubov transformation supposes the re- 
placement of the naive initial state by a transformed vacuum state, annihi- 
lated by a superposition of annihilation aa(k) and creation a^(— k) operators. 
For this purpose first we define the transformation 



(B.l) 



and then a new vacuum state |0 > 
a^(k)|0 >= 



70 



n 



aO 



aa(k) 



n 



130 



p„^a+(-k) 



10 >= 



(B.2) 



with p = C symmetric matrix. By means of the expectation value of 
aa{k)at(k') in the new vacua we define a matrix M^^ 



< 0|a„(k)a+(k')|0 >= {2rt)^6^{k - k')2y^n:^oM^^{k) 
which is related to the matrix p as follows 



M - pM^p^ 



I . 



(B.3) 



(B.4) 



Following the schema presented in Appendix A of [13] one sees that the 
cancellation of the dangerous terms we have mentioned above requires 

1 1 



Imp"^ 



130) 



rV:^(0) , 



(B.5) 
(B.6) 



So after p is known we have to determine the fluctuation integrals and 
the energy-momentum tensor in the Bogoliubov-transformed vacuum. But 
for the numerical implementation it is more convenient to redefine the mode 
functions 



(B.7) 



where the matrix is defined as N x N = M. Here f" are solutions of 
the same mode equations as are. For the fluctuation integrals and the 
energy-momentum tensor the expressions presented in Sec. 13.31 remain valid, 
with the mode functions /" replaced by the Bogoliubov transformed ones 
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C Perturbative Expansion 

The perturbative expansion of the mode functions for the case of coupled 
system in curved space-time follows a scheme introduced in Ref. |37|. We 
split the mode functions /"(r; k) into the free part and higher order terms 
represented by the reduced mode functions hflr; k) with the ansatz 

ftir;k) = e-^''-''-[fto + htir;k)]. (C.l) 

The mode equations fl3.25p are equivalent to the following integral equation 

/f (r; k) = e-^^-Vfo + f dr' A,,{r - r'; A;)V,,(r')4"(r'; k) , (C.2) 



with the potential Vij introduced in Eq. fl3.44p and with the retarded kernel 
of the free equation 

A,,(r - r'- k) = J2 ^©(^ - rO/fo/jlr'"^^^-^'^ " e-^''o(— ')] . (C.3) 
Then the reduced mode functions satisfy the integral equation 

htir- k) = r dr' J2 -7^/fo//o[e'^''"°''''''°^^"""'^ - e'(^"°-^^«)(^-^')] 

Jo ^ /30 

xyAr')[fko + Kir';k)] (C.4) 
and the corresponding differential equations become 

hf - 2zQ^ohr{r; k) = -Y, [v.^Z/o + (M^ - ml,6,,)h'^{T; k)] . (C.5) 

j 

The initial conditions ^M) and ^:Mf imply hf{0; k) = hf{0;k)=0. 

Now we expand the reduced mode functions /i" with respect to orders in 
potential Vij 

Kir) = hf\r) + hf'^ir) + hf\r) + . . . , (C.6) 

where /i"^"'' denotes the nth order in Vij. For our calculations we will need 
only the first and second terms of this expansion, therefore we introduce 

notation /i"^"''(r) for the sum over all orders in Vij starting with the nth 
order, i.e. /i"'^^'' = ^^^^ + /i"^^''. Using the iteration one can obtain reduced 
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mode functions of any order. In first order we have 



'JiOJ jOJ i 



P fa 
kO 



dT'Vjk{r')e 



Jo 



(C.7) 



The large-momentum behaviour of the Fourier-type integrals is analysed in 
Appendix To leading order in we find 



fiO^Pa 



T 



and 



Im h 



ail) 



[t; A;) ~ — / 
^ 2fi;30 Jo 



where we have introduced the notation 



(C. 



(C.9) 



(C.IO) 



Vq,/3 is symmetric in the indices, as is Vjk- 

For the fluctuation integrals fl3.36p we need to compute 

Re(/f/-) = f-f-+f-RehJ + f^,Reht + Reihfh 



T) 



(C.ll) 



For large momenta the first term is constant, of course. The second and third 
terms behave as see Eq. (IC.Sp . The term Re {h'^'"^^ h^^^"**) obviously 



contains Im/i" Tm/iJ which likewise behaves as see Eq. flC.Op . It 

would imply a nonlocal divergence. This contribution can be shown to cancel 
against a similar term in /j^Re /i"*'^'' + ffoR^ however. For the single 

field case this has been discussed in Ref. [37j (see there below Eq. (56)). 
The second and third term appear in the integrand of fluctuation integral 

as: 



E 



2n 



aO 



/foRe hf^ + /;oRe K 
Vpa 



EfarP 
JiOJ jO 

a,P 



a{l) 



2f^QO^/3o(^Q:0 + ^^0) 
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(C.12) 



where we have used Eq. fIC.Sp . 

We have aheady used the Wronski relation at r = when defining the 
initial conditions in Sec. 13.21 Now we will consider the Wronskian in gen- 
eral, which will be useful for simplifying the mode integrals occurring in the 
energy-momentum tensor 

W = frft - fr'ff = -^t^aoSa^ . (C.13) 

From this we obtain the relation 

^(^^.0 + n,o){ftoh^ + fl^K* + hTh^) = f^' - f^,hr' + K^h^' - hT'h^ . 

(C.14) 

For our purposes we will need only the a = [5 part of this relation, i.e.: 

n,o(2/rRe /i^ + = /°Im hf + Im h1*hf . (C.15) 

This relation can be used to simplify the kinetic term in the quantum part 
of the tt component of the energy-momentum tensor. Using the expansion 
in terms of reduced mode functions and Eq. (IC.lSp one finds 

E Cf^ = E ^ [^loUto - 2/roRe ht - K?) + \htf] . (C.16) 



2^00 2r2^_ 

a a 

This can be combined with 



E ^ = E ^ [f^' + ht + \ht\'] (C.17) 

a a 

to obtain 



When integrated over momentum the first term in the bracket, including the 
pref actor l/2flao leads to a quartic divergence. The second term is dominated 
at large momenta by the square of the imaginary part of /i"^^^ (seeEq. KJS)\\ ): 

up to terms of order This implies a logarithmic divergence. Therefore 
the subleading terms are not important, and we may replace this asymptotic 
estimate by 

2(r i^^'i' E AO o (o , Q 1 ^/^"^"^ • 

This has the advantage that it can be integrated in closed form and is of the 
same form as other logarithmically divergent contribution, see Eq. (IC.12p . 
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D Large-momentum behaviour of a Fourier 
transform 



When analysing the reduced mode function h°'^^\ Eq. (lC.7p . at large mo- 
mentum we encounter the expression 



(D.l) 



For a = (3 this integral has been analysed for large k in Ref. 
For a ^ (3 the analysis is somewhat more involved. The first exponential 
exp [i{Qao + ^/3o)(t — t')] will oscillate strongly when — )• oo. So this part 
of the integral will be dominated by the region around r = r', i.e., the upper 
end of the integration interval. The large-momentum behaviour can be anal- 
ysed by repeated integrations by parts, as for the case a = (3, see below. The 
second exponential exp [i{flao — ^/3o)(t ~ t')] will approach unity as — > oo. 
Indeed 

n^o - = - 7 (D.2) 

as /c — )• oo, SO the exponent goes to zero in this limit. 

Considering the two parts separately we find for the first one, via inte- 
grations by parts, 

1 







"aO + "00 



i 



K,.(r) - v;fc(0)e'(^-°+^^«)"] + . . . (D.3) 



and for the second one, by expanding the exponential 







2 



= / dr'V,,{r') - ^^^^^-^ / dr'V,,{r')ir - r') 

Jo "aO + "^0 Jo 

Z[\LaO + "^Oj Jo 

When contracted with fjofj^Q the second term on the right hand side does 
not contribute as Vai3 = /fo/fcoV?fc symmetric in a and (3. 
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E Dimensionally regulated integrals 

For the dimensional regularisation we need some identities [D 
d^-^k 1 



/30j 



-In 



167r2 
d^-^k 1 



+ 1 + 



/30 



In- 



/30 



167r2 



/i2 



/30 



+ 1 



aO 



TT 



327r2 



/i"' 2 



4-e): 



(E.l) 

(E.2) 
(E.3) 



where = | — 7 + ln47r. 

In the hmit m^Q — ?■ m^g the finite terms in the first of these equations 
reduce to 



hm 



/i2 



+ 1 + 



/30 



2 2 



In 



/30 



aO 



(E.4) 



F Counterterms 

Using the exphcit form of M-f^ the counterterm Lagrangian Eq. (14 .Sp can be 
decomposed into the usual mass, coupling constant and other counterterms 
as 

C^' = --{6m^ + 6^,R)a'<jy^-^6X,,<p^<f>'^ (F.l) 

Note that the first terms do not have a prefactor a", as it has been incorpo- 
rated into the rescaled fields (pi. 
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When we compare this to the expansion of the right hand side of Eq. (14. 5p 
we find 



5mf 



6Xi. 



6Xij 



6Z 
6a 
5A 



1 



167r2 
1 

16^' 
1 

1 



3A,,(6-l/6) + ^A,,(e,— 1/6) 



9A^ + E4 



^3\ ' 



i 



327r 



647r2 



4 



(F.2) 
(F.3) 

(F.4) 

(F.5) 
(F.6) 

(F.7) 

(F.8) 



These are in agreement with the counterterms found for the case of a single 
field We recall that the last three renormalisation constants are re- 

lated to renormalisation of the left-hand side of Einstein's equations, which 
have been redefined as renormalisation of the energy-momentum tensor, see 
Eqs. and f lMj) . 

The counterterms do not by themselves define the renormalisation. 
Choosing them as given above defines the MS scheme. They can be mo- 
dified by finite terms in order that the renormalised theory satisfies certain 
requirements. E.g., if we want our universe to become matter or radiation 
dominated at late times we will have to ensure that the effective cosmolo- 
gical constant remains zero. We will also want Newton's constant to retain 
its observed value. The details of such finite renormalisation depend on the 
specific model and usually involve tedious calculations. We do not address 
this issue here. 



G An instructive example 

The fluctuation integrals obviously play a central role in the one-loop appro- 
ximation to non-equilibrium dynamics. In particular it plays a role for fixing 
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the parameters in the equihbrium theory, one would not hke the renormah- 
sation conditions to be dependent on the expansion of the universe. As we 
have mentioned aheady the In a(r)- term in the finite parts is new here 
and suggests a logarithmic increase. It can become very large if the universe 
expands by many e-foldings. So it is interesting to analyse the behaviour of 
the subtracted integral J^"^ as well. We will do this for just one field, and 
setting for simplicity ^ = 1/6. Furthermore we will assume that the universe 
has ended up, at least approximately, in its ground state, i.e. that the clas- 
sical field (j) has reached its equilibrium value, which we assume to be and 
that the fiuctuations have almost died out, in a sense we will specify shortly. 
The mode equation for /(/c, r)is given by 

r + (A:' + -M'(r))/ = 0. (G.l) 

If the classical field has reached its minimum Jvi^ is still time dependent 

M'iT) ~ m^a' , (G.2) 

so the modes are still under the influence of a time-depending "potential" . 
At least for large k we can approximately solve the mode equation using the 
semiclassical approximation 

/(A;,t) ~C(k)(A;2 + M2(r))-V4exp(±i5) , (G.3) 

with 

S = I dr' ^ + M^{t') . (G.4) 

J rn 



So we have 



i/wr-^^^ffi=. (G.5) 

We will fix the constant C{k) by the assumption that the system is, at least 
as far as the high frequency modes are concerned, already in its adiabatic 
vacuum at tq. Then for sufficiently large k 



C\k)c^ ^k^ + M^{ro) . (G.6) 
Then with Ai^ = rn^tt^ and A^^(ro) = itlq we have 

V(t) = m^a^{T) - mO^ . (G.7) 
and the integrand of subtracted fluctuation integral becomes 
k^ 



47r2 



+ 



2 2 2 

ma — rriQ 



2 -L ml" 



(G.8) 
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Integrating from /c = to infinity we obtain 

with corrections due to the low-energy modes. The finite term J-"^* here is 
given by 

1 



167r2 



2 

2 2 1 2 

mam — 



(G.IO) 



where we have used the single-field version of Eq. (14.111) For the finite value 
T^"^ we therefore get 

^fin ^ ^ft ^ ^sub ^ ^ An _ 1 V (G.ll) 

The logarithmic dependence on a has disappeared. The factor is natural. 
In the equation of motion the fluctuation integral adds to 0^ = a^0^. 

As an application of interest we may consider the broken symmetry case. 
We have a Higgs field X; with a classical vacuum expectation value v. If 
its self-coupling is denoted by a and if we use the MS scheme the one-loop 
corrected equilibrium position is determined by 

< >= v^a^ + aJ^^"" . (G.12) 



With Eq. (IG.lip we see that the expectation value of x itself becomes inde- 
pendent of a, as it should. Were it not for the In a term in J-"^*, the vacuum 
expectation value would shift continuously with the expansion of the uni- 
verse. 

Of course there are corrections, but this is to be expected. In particular, 
if the fluctuations have thermalised, the vacuum expectation value of the 
Higgs field becomes a function of temperature. So our estimate would be 
strictly valid only after the temperature has reached 0. 

We would like to add that the heuristic estimates for the behaviour of the 
quantum fluctuations which we have used here are corroborated by numerical 
simulations. Actually they may be used for estimating the high momentum 
contribution to the subtracted fluctuation integral. 
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